Linear, diatomic crystal: single-electron states and 

large-radius excitons 

Vadym M. Adamyan* Oleksii A. Smyrnov'*' 
Department of Theoretical Physics, Odessa I. I. Mechnikov National University, 
2 Dvoryanskaya St., Odessa 65026, Ukraine 

January 2, 2009 

Abstract 

The large-radius exciton spectrum in a linear crystal with two atoms in the unit cell was obtained using 
the single-electron eigenfunctions and the band structure, which were found by the zero-range potentials 
(ZRPs) method. The ground-state exciton binding energies for the linear crystal in vacuum appeared to be 
larger than the corresponding energy gaps for any set of the crystal parameters. 

PACS number(s): 73.22.Dj, 73.22.Lp, 71.35. Cc 

1 Introduction 

The study of the quasione-dimensional semiconductors with the cylindrical symmetry became an urgent problem 
as soon as investigations of semiconducting nanotubes had been launched. One of the most important trends 
of research in this field is the study of optical spectra of such systems, which should include the exciton 
contributions [I]- [9]. Evidently, the quasione-dimensional large-radius exciton problem can be reduced to the 
ID system of two quasi-particles with the potential having Coulomb attraction tail. Due to the parity of the 
interaction potential the exciton states should split into the odd and even series. In [TO] we show that for the 
bare and screened Coulomb interaction potentials the binding energy of even excitons in the ground state well 
exceeds the energy gap (in the same work we also discuss the factors, which prevent the collapse of single- 
electron states in isolated semiconducting single- walled carbon nanotubes (SWCNTs). But the electron- hole 
(e-h) interaction potential and so the corresponding exciton binding energies may noticeably depend on the 
electron and hole charge distributions. So it is worth to ascertain whether the effect of seeming instability of 
single-electron states near the gap is inherent to the all quasione-dimensional semiconductors in vacuum or it 
maybe takes place only in SWCNTs for the specific localization of electrons (holes) at their surface and weak 
screening by the bound electrons. That is why we consider here the simplest model of the quasione-dimensional 
semiconductor with the cylindrical symmetry, namely the linear crystal with two atoms in the unit cell. The 
electrons (holes) in this crystal are simply localized at its axis. 

The aim of this work is only a qualitative analysis of the mentioned effect. For study of electron structure 
of concerned ID crystal we apply here the zero-range potentials (ZRPs) method [II],[T2] (see section 2). The 
matter is that results on the band structure and single-electron states, obtained by this method for SWCNTs 
in |13|.|14|. appeared to be in good accordance with the experimental data and results of ab initio calculations 
related to the band states. For certainty we use the linear crystal parameters (the electron bare mass, lattice 
parameters) taken from works on nanotubes [13j.[14j. In section 3 we obtain the e-h bare interaction potential 
and that screened by the crystal band electrons, and then the large-radius exciton spectrum for the linear crystal 
in vacuum. All these data are used in section 4, where we present results of calculations for the crystal with 
different lattice periods (it also means different band structures). As it turns out, the binding energy of even 
excitons in the ground state well exceeds (~ 2 — 5 times) the energy gap for the linear crystal in vacuum and the 
screening by the crystal band electrons is negligible. Note, that this result was obtained within the framework 
of exactly solvable ZRPs model with feasible parameters. Therefore, the mentioned instability effect may take 
place not only for the considered simplest case, but, most Hkely, also for other quasione-dimensional isolated 
semiconductors in vacuum. 

*E-mail: vadamyan@onu.edu.ua 
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2 Single-electron band structure and eigenf unctions of band elec- 
trons 



We have obtained the single-electron states in the linear crystal using the zero-range potentials (ZRPs) 
method The main point of this method is that the interaction of an electron with atoms or ions 

of a lattice is described instead of some periodic potential V{r) by the sum of Fermi pseudo-potentials [11] 
(1/a)^; d{pi){d/dpi)pi (pi = \i — r;|, r/ are the points of atoms location, a is a certain fitting parameter) or 
equivalently by the set of boundary conditions imposed on the single-electron wave function ip at points r; : 

Pi^o [dpi J 

The electron wave functions satisfy at that the Schrodinger equation for a free particle for r 7^ r;. Therefore 
we seek them for the linear crystal in the form: 

i^[Pn,Pn)^ 2^ An + 2^ "b > (2-1) 

n— — oo ^ Ti— — C30 ^ 

where indices A and B denote two monatomic sublattices of the diatomic lattice, = | ^ '"n | ^'^'^ Pn ^ 
|r — rj^|, n numbers all the sublattices points, k = ^j2m\,\E\/h, i? < is the electron energy and mb is the 
bare mass. For certainty, following [13], [2] we take from now on mb — 0.415mc and the ZRP parameter 
a = ^2mb|i?ion|/^, where E^ion is the ionization energy of an isolated carbon atom. By [13j.|14|. with these a 
and mb ZRPs method reproduces single-electron spectra of such quasione-dimcnsional structures as SWCNTs 
within an accuracy of existing experiments. One can take infinite limits for the series in (|2.ip even for the finite 
crystal, because terms of these series decrease exponentially with increasing of n. 

According to the ZRPs method the wave functions (|2.ip should satisfy the following boundary conditions at 
the all sublattices points: 

hm (A (pjV) (r) -f a (pjV) (r)l = 0, (2.2) 

here i = {A, B} according to each sublattice. 

Further we suppose that the linear crystal lies along the z-axis, thus — nde^ and rj^ — [nd + a)ez, 
where is the z-axis unit vector, a is the distance between atoms in the unit cell of the crystal and d > 2a 
is the distance between the neighbour atoms in each sublattice. Note, that d = 2a corresponds to the metallic 
monatomic crystal and for the case d < 2a the smallest distance between atoms in the crystal is d ~ a < a. 

Substituting (|2.ip to (|2.2p and applying the Bloch theorem (A„ = Aexp{iqdn), Bn = Bexp{iqdn), q is the 
electron quasi-momentum) we get two equations for amplitudes A, B: 



AQ1+BQ2 = 0, 
AQl + BQi = 0, 



(2.3) 



where ^ 

Qi{k, q) — a — —\n {2[coshKd — cosqd]) , (2-4) 

X exp{-K\nd + a\+iqnd) 

^ ■ ^ ^ ^ 



n— — 00 



Setting d = ja: 

1 

Iff exp[— Kfl] a;''~^exp[Ka 



'^^ = a i U--^exp[d(i9-.)] + expKig + .)]-.. ) ^'-^^ 


for each real j > 2. 

From (|2.3p we get two equations, which define the band structure of the crystal: 

Qi(Ki,g)-|Q2(Ki,<7)| =0, (2.7) 

Ql(K2,g) + |Q2(K2,(7)| =0. (2.8) 
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Equation (|2.7p defines the conduction band and equation l|2.8p defines the valence band (see section 4, figure 1). 
So the electron and hole effective masses can be simply obtained from (|2.7p and (|2.8p . respectively. 

Further, using the Hilbert identity for Green's function of the 3D Helmholtz equation, we obtain the nor- 
malized wave functions (j2.ip : 

A{K,q) 1 exp{— K\r — nde^l + iqnd) 

^-'^^^ =^r [ ^ 

^ (2-9) 

Qi s—^ exp(— K|r — (nd + aje^l + iqnd) \ 
n~"oo \r~{nd + a)e,\ J ' 

where L is the crystal length and A{k, q) is the normalization factor: 

\ 1/2 

1 f Kd cosh Kd — cos qd ^ 



^ ''^'^ 2 \ TT sinh Kd — 5Ry 

and 



^ (exp[— igd] sinh ku + sinh K[d — a]) . 



3 Exciton spectrum and eigenfunctions. Bare and screened e-h in- 
teraction 

Using the same arguments as in the 3D case one can show (see, for example lOj), that the wave equation 
for the Fourier transform (j) of envelope function in the wave packet from products of the electron and hole 
Bloch functions, which represents a two-particle state of large-radius rest exciton in a (quasi)one-dimcnsional 
semiconductor with period d, is reduced to the following ID Schrodinger equation: 



2fi 



Eg, —oo < z < oo, 

where n is the e-h reduced effective mass and V{z) is the e-h interaction potential: 

V{z) 



(3.1) 



//((xi 



{{Xl - + (2/1 - 2/2)2 + {z + Z,- Z2)2)l/2 

X |uc;K,Vd(»^i)n"v;K.^/d(^'2)pdridr2, 
= E2 X (0 < z < d). 

Here Uc,v;K,g('") are the Bloch amplitudes of the Bloch wave functions i/'c,v;K.g('") = exp(igz)Mc,v:K.g('") of the 
conduction and valence band electrons of the linear crystal, respectively. Using the actual localization of the 
Bloch amplitudes at the crystal axis, after several Fourier transformations and simplifications we adduce the 
e-h interaction potential to the following form: 



^1 ,r2 (^) — 



Ji(fc)Ji(fcr2/ri) f k 



k^ 



-{\d- z\ + \d+ z\-2\z\ 



(3.2) 



- exp 



-\d-z\ 



exp 



\d + z\ 



2 exp 



dfc, 



where J is the Bessel function of the first kind and ri (r2) is the radius of the electron (hole) wave functions 
transverse localization 

1/2 



2 J "^20 y l"c,v:K,g(2,'r2D)|^dzdr2D 1 
^ E2 / 



where r2D is the transverse component of the radius-vector, q — n/d and k — ki^2 {T^/d) correspond to the 
conduction and valence bands edges at the energy gap (according to (|2.7p and (|2.8p . respectively). Equation (|3.ip 
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with the potential given by (|3.2p defines the spectrum of large-radius exciton in the linear, diatomic crystal if 
the screening effect by the crystal electrons is ignored. Actually, the screening of the potential (|3.2p by the band 
electrons is insignificant. 

Indeed, following the Lindhard method (so-called RPA), to obtain the e-h interaction potential (p{r), screened 
by the electrons of linear lattice, let us consider the Poisson equation: 

-A<^(r) = 4^(p-'(r)+p-^(r)), (3.3) 

where r is the radius- vector, p°^^{r) is the density of extraneous charge and /9'"'*(r) is the charge density induced 
by the extraneous charge. 

By p.3|) the screened e-h interaction potential may be written as: 



^{r) J {p'^'^'ir') + p'"'*(r')) G(r, r')dr', (3.4) 

E3 

where G{r, r') — l/(47r|r — r'\) is Green's function of the 3D Poisson equation. 

Let E'^{q) and -0° ^(r) = exp(igz)u° ^(r) be the band energies and corresponding Bloch wave functions of 
the crystal electrons and E{q), ipK,q{r) be those in the presence of the extraneous charge. Then 

r'^ir) ^ -e^ [f{Eiqm,Jr)\^ - /(i?°(g))K,(r)p] , (3.5) 
q 

where / is the Fermi-Dirac function. Using the transverse localization of the Bloch wave functions near the 
crystal axis, we get in the linear in ip approximation: 

L 

P"^'^{z',r'2B) ^ -e^y2~^ / r^^-,'^2,q'i^',r2i))u*.^^ {z,r2i))dr2i)(p{z)exp[iz{q' - q)]dz 

„ „, ^g;q,q' J J (3.6) 
9,9 E2 

X <;k2.<j'(^'' r2B)uc:Ki,q{^', t^d) exp[iz'(g - q')], 

where Eg.q,qi — Ec{q) — E^{q'). Here and further (f{z) is the e-h interaction potential averaged in E2 over the 
region of the Bloch wave functions transverse localization and over the lattice period d along the crystal axis. 
Due to the periodicity of the Bloch amplitudes p™'* may be written as: 

ind/ ^ C{q,q';d) 

L ^ £^g;9,9' (3.7) 
X K:K2^q'{r')u^,^,,,q{r') exp[iz'(q - q')], 

where 

d 

C{q,q';d) = J J u*.^^ g{z,r2D)uy-n2,q'{z,r2D)<ir2Ddz 
E2 

and N is the number of unit cells in the crystal. 

Further, after several transformations we obtain from p.4p and (j3.7p the one-dimensional Fourier transform 
of the potential ip: 



= 7(^7' 

^ , e^N^ ^f' \Ciq,q-k-d)\\ ~ 2sin(fcrf/2) ^^'^^ 

where (po is the Fourier transform of the averaged electrostatic potential induced by p°^^ and Ko{k) is the 
modified Bessel function of the second kind averaged over r2D and r'^j-, in the region of the Bloch wave functions 
transverse localization in E2 , namely 

Ko{k) = , ^ / / Koi\k\\r2D - r'2j,\)dr2Bdr'2^, 
['Krir2Y J J 

E2' = (0 < r2D < n) X (0 < /3 < 2tt). 
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In the long- wave limit we get: 



2 l2 



\Ciq,q-k;d)\t^,^\U{q;d)\'k 

d 



(3.9) 



Ea 

Using of the Schrodinger equation for the orthogonal Bloch wave functions yields 

d 



U{q;d) = 



//c...(..r.o)|*....(..r.o)d.,od. 

Eo 



(3.10) 



Hence, in the long-wave limit the screened quasione-dimensional electrostatic potential induced by a charge 
eo, distributed with the density: 



eo 



nR'^d 
0, X < xo 
X > xo 



(6 [z + d/2] -e[z- d/2]){Q [r2D] - 6 - R])-. R > Q, 



&[x-xo] = I J' 
in accordance with (|3.8p and p.lOp , is given by the expression 



ith 



9d 



TTd^ 



Seori f {l/k^)sin^{kd/2ri)Ko{k/ri)cos{kz/ri) 



1 + fjdikn/d) sm{kd/2ri)Ko{k/ri) 



dk 



V Trrimb 



Tv/d 



-TT/d 



V'c;ki ,q 



d_ 
'd'z 



''Pv;K2,q 



dq. 



According to equation (|2.9p the dimensionless screening parameter may be also written as: 



gd 



( 16e 



\ hdr\ 



2 ^ Ai 



Al{Kuq)Al{K^,q) 



{^i{q)^^l{q)f 

^ ^ Ql(Kl,g)Ql(^2,g) 

Ql{K.i,q)Q2{K,2,q) 

K2 

l(^l,g) [ X, 

^7 



Qi(k, g)dK 



Ql(«2,g) 
Q2{K2,q) 



Q2('«, 9)dK 



dg. 



(3.11) 



(3.12) 



(3.13) 



Note, that ki and K2 are the implicit functions of q defined by (|2.7p and (|2.8p . respectively. 

It appears, that calculated according to (|3.13p for d varying in the interval [2.1a, 3a] are about 10" 



4 Discussion. Stabilization of single-electron states 

Using equations (|2.7p and (|2.8p we obtained the band structure (see figure 1) and the electrons and holes effective 
masses for the linear crystal of dimers for different values of the ratio j — d/a of its period d and the distance 
a between atoms in dimers. Besides, using wave equation p.ip and potentials p.2p and (|3.1ip we found the 
large-radius exciton energy spectrum in the crystal for the bare e-h interaction and e-h interaction screened by 
the bound electrons of the crystal. We present here results for the crystal with j G [2.1,3]. Contrary to the 
single-band metallic crystal with j = 2, the crystals with j > 2 are semiconductors with band gaps varying from 
zero to the difference between the electron levels in an isolated dimer. Particularly, the crystals with j = 2.001 
and realistic values of a and a (as in nanotubes and some ID polymer chains) are narrow-gap semiconductors, 
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Figure 1: The band structure of the linear crystal with parameters: j — 2.1 (dashed line), j — 2.5 (dot-dashed 
line) and j = 3 (solid line); q in units of n/d. The upper and lower bands correspond to equation (|2.7p and p.8|) . 
respectively. 



Table 1: Band gaps Eg and reduced effective masses /i according to ()2.7|) . (|2.8|) : radii of the electrons and 
holes transverse localization ri and r2, respectively; screening parameter according to ()3.13p and the exciton 
binding energies £ of the even and odd series for the linear, diatomic crystal according to equation p.ip with 
potential (|3.2p for different values of the ratio j — d/a. 



j 


Eg (eV) 


/i (TOo) 


ri (nm) 


r2 (nm) 


9d (10-^) 


iSoieven (eV) 


i^^i;odd (eV) 


^0;cven/ -^g 


2.1 


1.4422 


0.041 


0.070 


0.0611 


0.7235 


-6.90 


-0.5939 


4.7845 


2.3 


3.3146 


0.125 


0.080 


0.0569 


2.4716 


-8.9992 


-2.0631 


2.715 


2.5 


4.403 


0.2199 


0.088 


0.0549 


2.7036 


-9.5352 


-3.4812 


2.1656 


3 


5.6281 


0.5665 


0.0994 


0.0551 


1.1206 


-9.6588 


-5.8421 


1.7162 



in which excitons may possess binding energies about ^ 10 meV, but the crystals with j = 2.1 are already 
wide-gap (^ 1 eV) semiconductors with strongly bound e-h pairs, and the crystals with j ~ i are almost flat 
band semiconductors, but their electrons and holes at the energy gaps {q = tt/c?) still have the finite effective 
masses (these electrons and holes form the excitons in the crystals). For certainty, the distance a we have 
chosen equal to the graphite in-plane parameter 0.142 nm. The ZRP interaction parameter a = 11.01 nm^^ 
corresponds to the ionization energy of an isolated carbon atom (i^ion = 11.255 eV). 

As one can see from table 1 the obtained from (j3.13p dimensionless screening parameter gd 1 for the all 
considered values of j. So, it turns out that the screening of the e-h interaction potential by the band electrons 
in the linear, diatomic crystal may be ignored. This result could be expected since the considered model of 
linear crystal is close to that of the electron gas confined to a cylindrical well. In the latter case, for which the 
separation of the angular variables takes place, the states with different quantum numbers m of the angular 
momentum play the role of electron bands. Accordingly, the matrix element | {ipc \d/dz \ ipv) \ from (|3.12p for the 
direct transitions between bands with different m appears to be identically equal to zero. This is why only the 
binding energies of excitons with unscreened interaction potential are listed in table 1. 

To obtain estimates of the main linear crystal characteristics we considered several limiting cases. In the 
case of d ^ a (or j 3> 1) and a = const equations (|2.7p and (|2.8p can be reduced to a — Ki_2Texp[— Ki^2a]/a = 0, 
thus bands become flat (ki.2 do not depend on q) and the band gap tends to the finite value (fi^/2mb)(K2 — ) 
(for a ~ 0.142 nm it is about 6.3 eV), hence the reduced effective mass and exciton binding energy tend to 
infinity, while the exciton radius Toxc ~ h^/fj,e^ tends to zero. Therefore, the large-radius exciton theory is 
actually appropriate (^oxc S> a) for excitons in the linear, diatomic crystal only when its period d runs the 
interval (2a, 2.4a) (e.g., Toxc is ~ 9a for j = 2.1 and ^ 2a for j — 2.4). If d = const, but a — > 0, the conduction 
band moves to the region of positive energies and at some critical value of a disappears, while the valence band 
shifts correspondingly to the region of deep negative energies. 

Table 1 shows that the ground-state exciton binding energies for the linear, diatomic crystals with any value 
of the ratio j are greater than the corresponding energy gaps. Note, that according to the same calculations, 
but with the bare mass mb = rric, the ground-state exciton binding energy for the linear crystal in vacuum 
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appears significantly greater than the energy gap. We should note also, that the ground-state binding energies 
of excitons in the linear crystal with different periods d in vacuum, calculated using the ID analogue of the 
Ohno potential |15j instead of potential (j3.2p . remain greater than the corresponding energy gaps. Particularly, 
for d — 2.3a calculations with the ID unscreened Ohno potential with the energy parameter U taken from [TB] 
{U ~ 11.3 eV) and [17] {U = 16 eV) give the ground-state exciton binding energies £o;cvon = 5.90 eV and 
^^O:ovcn = 7.63 eV, respectively, while Eg = 3.31 eV (see table 1). Thus, all calculations made on the base of 
solvable zero-range potentials model indicate the instability of the single-electron states in the vicinity of the 
energy gap with respect to the formation of excitons. This might be a shortage of this model, but it is worth 
mentioning that results obtained on one-particle states in real ID systems, like SWCNTs [13], [13], within its 
framework agree with existing experimental data in limits of accuracy of the latter. 

The stability of single-electron states of ID semiconductors with respect to the exciton formation in vacuum 
can be explained by bringing multi-particle contributions into the picture. With the advent of some number of 
excitons in the quasione-dimensional crystal the additional screening appears, which is caused by a rather great 
polarizability of excitons in the longitudinal electric field. This collective contribution of born excitons into the 
crystal permittivity returns the lowest exciton binding energy £o:ovcn into the energy gap and so blocks further 
spontaneous transitions to the exciton states. To show this let us consider the model of linear crystal immersed 
into the gas of excitons with dielectric constant Soxc confined to the region of linear crystal carriers transverse 
localization, namely: cylinder with radius ri and axis coinciding with that of linear crystal (from now on, for 
estimates, we assume that electron and hole have the same transverse localization radius). In this case it is easy 
to show that the e-h interaction potential is given by: 



IGe^ri fsm^{kd/2ri)cos{kz/ri) 



rd^ 



.k^ 



X 1 



2K^{k)h{k) 



(4.1) 



k{e,^,Ko{k)h{k) + Ki{k)h{k)) 



d/fc, 



where /; and Ki are the modified Bessel functions of the order i of the first and second kind, respectively. 
Further, like in |18j , we use the known elementary relation between the permittivity of exciton gas and its 
polarizability a in the direction of linear crystal 



1 + Ana. a — 2e n 



E 



(^o|r|vl/fc)P 
£() — £k 



where n is the bulk concentration of excitons, ^'o and £q are the exciton eigcnfunction and binding energy, 
which correspond to the ground state, and and are those, which correspond to the all excited states of 
exciton. Thus, the upper and lower limits for a are: 



2^'" -|(*o|r|*i)P <a< |(*o|r|vI/, 



£o — £i 



2e^n 
So — £i 



where and £i correspond to the lowest excited exciton state. Hence, one can obtain the upper and lower 
limits for n: 



1 £o;e 



£l;o. 



dd 



4-77 



2e2 



< n 



< 



^cxc 1 ^0;cvcn £l\0' 



dd 



477 



2e2 



z(j)o{z)(l)i{z)dz 



(4.2) 



where each (j) is the component of Fourier transform of the corresponding exciton envelope function, it depends 
only on the distance z between the electron and hole. At that, (j)o is the even solution of (|3.ip with potential (|4.ip . 
which corresponds to the exciton ground state and satisfies the boundary condition (j)'{0) = 0, and is the odd 
solution of the same equation, which corresponds to the lowest excited exciton state and satisfies the boundary 
condition 0(0) = 0. 

Varying Scxc in (|4.1[) substituted into the wave equation (j3.1[) one can match £o;ovcn to the forbidden band 
width. Further, f i;odd can be obtained from the same equation with the fixed Sexc and with the corresponding 



7 



boundary condition. All these magnitudes allow to calculate from (|4.2p the rough upper and lower limits for 
the critical concentration of born excitons ric, which is sufficient to return £o:cvon into the energy gap. Further, 
knowing ric we can calculate the shift of the forbidden band edges, which move apart due to the transformation 
of some single-electron states into excitons. This results in the enhancement of energy gap 



{nhnc)"^ 
2n 



SE, (4.3) 



like in [19j and [20] . Here rif. — rif-Trrf is the linear critical concentration of excitons, and ri is the radius of 
electron wave functions transverse localization at the linear crystal axis. 

In accordance with (|4.2p the described model yields ric ^ 180 fim^^ (~ 3% of the atoms number in the 
crystal) and ~ 400 /im~^ {^-^ 7%) for the linear crystal with j = 2.1 and j = 2.3, respectively, while by (|4.3p the 
corresponding SEg are ~ 300 meV and ~ 500 meV in the same order. Here, however, we should mention that 
for SWCNTs both the measured in [19], [20] and estimated in the same manner [18] values of SEg/ Eg appeared 
to be two-four times less. 

Note, finally, that the instability of the single-electron states weakens or disappears for linear crystals 
immersed into dielectric media. As it was shown in [9J by the example of the poly-para-phenylenevinylene ID 
chain or in [K] , [17] , [21] by the example of SWCNTs the environmental effect may substantially decrease the 
excitons binding energies. Indeed, for the linear crystal in media even with permittivity about e ^ 2 — 3 (e.g., 
like in [IB] or 07]) the ground-state exciton binding energy becomes smaller than the energy gap. 
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